Isospin-projected nuclear level densities by the shell model Monte Carlo method 
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We have developed an efficient isospin projection method in the shell model Monte Carlo ap- 
proach for isospin-conserving Hamiltonians. For isoscalar observables this projection method has 
the advantage of being exact sample by sample. The isospin projection method allows us to take 
into account the proper isospin dependence of the nuclear interaction, thus avoiding a sign problem 
that such an interaction introduces in unprojected calculations. We apply our method in the cal- 
culation of the isospin dependence of level densities in the complete pf -\- gg /2 shell. We find that 
isospin-dependent corrections to the total level density are particularly important for N ^ Z nuclei. 



The level density is among the most important statisti- 
cal properties of nuclei. It is required for the calculation 
of transition rates through Fermi's golden rule and of nu- 
clear reaction rates through the Hauser-Feshbach theory. 
Applications in astrophysics include estimates of neutron 
and proton radiative capture rates in the s, r, and r-p pro- 
cesses [H, 0- 

Most theoretical approaches to the level density are 
based on the Fermi gas model, i.e., Bethe's formula 
This formula describes the level density of a non- 
interacting many-fermion system in terms of the single- 
particle level density parameter a. Shell corrections 
and two-body correlations (e.g., pairing correlations) are 
taken into account empirically by shifting the ground- 
state energy by an amount A. The resulting backshifted 
Bethe's formula (BBF) for the total level density at ex- 
citation energy i?^; is given by 
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where t represents the Lang-LeCouteur modification 
defined \>y at^ - t ^ - l\. The BBF describes well 
the experimental level densities of many nuclei if the pa- 
rameters a and A are fitted individually for each nu- 
cleus 0. Although their global systematics have been 
studied extensively, these fitted parameters exhibit a sig- 
nificant dependence on the nucleus under consideration. 
Since nuclear level densities increase very rapidly with 
i?^:, an accurate determination of a and A is crucial for 
reliable calculations of level densities and reaction rates. 

As the basic microscopic model of nuclear structure, 
the interacting shell model takes into account both shell 
effects and correlations. This model has been success- 
ful in describing the low-lying states of many nuclei and 
is also an attractive framework to study microscopically 
thermal and statistical properties of nuclei. However, the 
large dimensionality of the many-particle model space 
has limited the applicability of direct diagonalization 
methods of the shell model Hamiltonian in medium- and 
heavy-mass nuclei. This limitation has been overcome us- 
ing the sheU model Monte Carlo (SMMC) approach 0,0]. 
While the SMMC method cannot provide a detailed spec- 



trum of the many-particle Hamiltonian, it enables cal- 
culation of thermal and statistical properties in model 
spaces that are many orders of magnitude larger than 
those that can be treated by conventional methods. 

We have developed a method to calculate level den- 
sities in the SMMC approach, and applied it to nuclei 
in the mass range A 50 — 70. We found good agree- 
ment with experimental data without any adjustable pa- 
rameters ■ Fermionic Monte Carlo methods often 
suffer from the so-called sign problem, which leads to 
large statistical errors and a breakdown of the method 
at low temperatures. However, the dominating collec- 
tive components of the nuclear interaction [ll| usually 
have a good sign, and general effective interactions con- 
taining small terms that have a bad sign can be treated 
with the method introduced in Ref. [7|. In the calcula- 
tion of level densities the sign problem was circumvented 
by constructing a good-sign interaction in the p/ -I- 59 /2 
shell [1| that correctly includes the dominating collective 
components of realistic effective interactions. 

In the SMMC approach, thermal observables are aver- 
aged over all possible values of the good quantum num- 
bers. However, many of the applications require the de- 
pendence of thermal observables on the good quantum 
numbers. This dependence is determined in the SMMC 
approach by projection methods. Parity projection was 
implemented in the study of the parity dependence of 
level densities [I,@,[i3,[l3|, and a spin projection method 
was introduced to determine the spin distribution of nu- 
clear levels |14j|. The particle-number reprojection tech- 
nique has been developed [l3] to facilitate the systematic 
studies of level densities for a large number of nuclei. 

The isospin T is approximately a conserved quantum 
number in nuclei. Reliable treatment of the isospin de- 
pendence of the level densities could be important in 
N ^ Z nuclei because levels with different T values lie 
close in energy. In our SMMC studies of level densi- 
ties we have used a good-sign Hamiltonian to keep the 
statistical errors small. While this Hamiltonian contains 
the dominating collective components of nuclear effective 
interactions, it might not describe properly the energy 
differences of states with different isospin values. Such 
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discrepancies can be corrected by including an isospin- 
dependent interaction 15, l^l- However such an interac- 
tion has a bad Monte Carlo sign and was treated pertur- 
batively in Refs. [H, [l6| . Here we overcome this problem 
by introducing an exact isospin projection method in the 
SMMC approach. This T-projection enables us to de- 
termine the isospin dependence of nuclear levels and to 
include accurately an isospin-dependent interaction term. 
We apply the method to the calculation of level densities 
in the complete pf + Ogg /2 shell. 

The SMMC method is based on the Hubbard- 
Stratonovich (HS) representation of the imaginary-time 
propagator e~^^ describing the Gibbs ensemble of a nu- 
cleus with Hamiltonian H at inverse temperature /3. In 
the HS decomposition e^'^^ — J D\cr]GaUa, where is 
„„„ - " " d.Th„(T) ^j^^ prop- 
agator of the one-body Hamiltonian h^^ describing non- 
interacting nucleons moving in time-dependent auxiliary 
fields cr(T) (T denotes time ordering). The canonical ex- 
pectation value of an observable O is then given by 



(O) 
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where Tr^i^T^^ denotes a trace in the subspace of the 
fixed mass number A and isospin component T^. Since 
A = N + Z and = (A^ - Z)/2, this trace is equiv- 
alent to a trace at fixed neutron number N and proton 
number Z. We have also used the notation {Xcr)w = 
J D[a]W{a)X^/ J D[cr]W{a) with W{a) = GaTTA,T,U^- 
Here and in the following we assume TTA,T^Ucr > 
(i.e., we assume a good sign interaction and an even- 
even or = Z nucleus). Otherwise the sign function 
'^^A,T;,Ua/\TTA,T^Ua\ has to be included. In the Monte 
Carlo method we choose M samples of the a fields accord- 
ing to the weight function W{a) (each sample denoted by 
CTfc), and estimate {Xcr)w « J2k^'^k/^- 

Let us consider a model space for a set of isobars (i.e., 
nuclei with fixed A but different values of T^). Since the 
isospin forms an su{2) algebra, isospin projection gen- 
erally requires a three-dimensional integration over the 
Euler angles in isospin space (in analogy with angular 
momentum projection). However, this integration is time 
consuming and in the following we describe a more effi- 
cient method. We assume the nuclear Hamiltonian H to 
be isospin invariant (i.e., [H, T] =0), in which case both 
T and are good quantum numbers and the isospin 
multiplets (at fixed T) are degenerate with respect to T^. 
The subspace with fixed (i.e., fixed nucleus) contains 
all energy eigenstates with T > \Tz\. In the following we 
assume N > Z without loss of generality, and denote by 
To = {N — Z)/2 > the value of Tz for a specific nucleus. 



If the operator X is an isoscalar, we can decompose its 
trace for the specific nucleus as 



Tr^,' 



T>Tc 



TtatX 
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where Tr^i^^ represents the trace for fixed A and isospin 
T (not including degeneracy with respect to T^). The 
trace of X for a specific isospin T = Tq is then obtained 
by 



Tta,t=ToX = Tta,t^=ToX — Tta,t^=To+iX ■ 



(4) 



We note that Tz = Tq + 1 corresponds to a nucleus with 
(iV-|-l) neutrons and (Z — l) protons. In the level density 
calculations we use X — He~^^ which is an isoscalar. 
If X is not an isoscalar, a T^-dependent factor is required 
in Eq. (O and hence in Eq. ([4]). 

An isospin-invariant two-body interaction can always 
be decomposed into a sum of squares of isoscalar one- 
body operators 0]. Therefore , the one-body Hamilto- 
nian ha and the propagator Ua (in the corresponding HS 
representation) remains isospin invariant for any config- 
uration of the (7 fields. Equation ([4]) then holds with 
X = OUcr as long as O is an isoscalar observable 

TTA,T=ToiOUa) = TTA.T^=T„iOUa)~TTA,T^=To + l{OUa) ■ 

(5) 

Thus, T projection can be carried out for each sample by 
using the simpler Tz projection. Equation ((S]) guaran- 
tees that T projection is implemented exactly sample by 
sample. An equation analogous to Eq. (j4|) was used for 
angular momentum projection with the scalar operator 
X — Oe~^^ However, in contrast to the isospin 

projection, in the HS transformation is not a scalar 
under spatial rotations so the equation analogous to ([5|) 
does not hold sample by sample. We therefore expect 
that the present T projection method leads to statisti- 
cal errors that are typically smaller than the statistical 
errors in the angular momentum projection of Ref. 14 1 . 

Based on Eq. (O, isospin projection can be im- 
plemented using the particle-number reprojection tech- 
nique [10], in which the Monte Carlo sampling is done 
for a reference nucleus {N, Z) and thermal observables 
are evaluated by reprojection on a nucleus {N',Z'). In 
particular, the ratio between the T- and T^-projected 
partition functions is given by 



Za,T=To (/?) _ ^ _ / TT^A,T,=To + lUa 
Za.T,=To{P) \ TrA,T,=To^^o- 



(6) 
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Similarly, the expectation value of an observable O for 
isospin T = Tq can be calculated from 
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{0)a,t^ 



T^A,T,=To{OU^) ' 



/ TrA,T,=To + iiOU^) 



T^A,T2=TQ + lUcr 

Tta,t^=t„U^ 



1 - 



TrA,T^=To + lU^ 

TlA.T,=ToU„ 
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The Monte Carlo sampHng in Eqs. ([6]) and ([7]) is carried nucleus A,To + 1) is then used to calculate {0)a,t=To- 
out according to the weight function W{a) for the nucleus Isospin projection on T = Tq + 1 can be implemented 
A,To. Equation ([7]) (that includes reprojection on the within the same Monte Carlo sampling using 



A,T=To + l 



TlA,T,=To + l(OU^) TrA.T,=To + l(U^) ' 
TrA,T,=To + l{U^) ' TvA,T,=ToiU„) /yy 



Tl-A.T,=To + 2iOU<r) 
TrA.T,=To+2{U^) 



Tr 



A,T^ = T„ + 2 



Tr^,T, = To((7„) 



w 
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TrA,T,=To(C/^) 



A,T^=To + l 



/ TrA,T,=To+2(Ua) 
\ Tia.T^^TqU^ 



(8) 
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The projection on higher values of T is carried out in 
a similar manner. Thus we can project on all possible 
isospin values T using the same Monte Carlo sampling. 
In particular, we can calculate the level densities of iso- 
bars with a single Monte Carlo sampling process as in 
the number reprojection method of Ref. [lO|. 

The good-sign Hamiltonian of Ref. [1] , while reproduc- 
ing the proper collective features at fixed isospin T, does 
not necessarily reproduce the correct isospin dependence 
of energy levels. This isospin dependence can be particu- 
larly important for N ^ Z (i.e., ~ 0) nuclei, in which 
number of levels with T ^ Tq = {N-Z)/2 and T = To + l 
are comparable even close to the ground state. A simple 
way to account for the proper isospin dependence of the 
nuclear interaction is to add an appropriate function of 
T^, /(T^), to the effective Hamiltonian. The simplest 
such function is /(T^) — aT^ (a is a constant), as in 
the modified surface delta interaction 17| . This term 
is repulsive {a > 0) and leads to a sign problem in the 
SMMC method. In Refs. [H, such a T-dependent 
interaction term was treated perturbatively. However, 
the criterion for applicability of perturbation theory does 
not always hold. As we explain below, the T projection 
method enables us to account exactly for an arbitrary 
function /(T^) in the Hamiltonian. 

To determine the proper isospin dependence in the 
Hamiltonian, we extract from experimental data the ex- 
citation energy (T) of the lowest level for each isospin 
value T. When not directly measured, the experimen- 
tal values of Ex{T) can be obtained from the measured 
masses of the corresponding isobars, together with esti- 
mates of the Coulomb energy differences among them. 
For the latter we use Ec{A, Z)-Ec{A,Z -1) = |(2Z- 
l)e^/Rc with Rc = l.2AA^/'^ im.. We then shift the cal- 
culated value of E{T) (the lowest energy of isospin T), 



by 5E{T) - E{Ta) + E,[T) - E{T) (for T > To), so 
as to reproduce the experimental value of Ex{T). This 
shift defines the function /(T^) in the effective Hamilto- 
nian. The modified T-projected level densities are then 
determined by shifting the excitation energy at each T 
by SEiT). 

We first apply the isospin projection method to the 
nucleus ^^Cu. The effective interaction consists of 
T — 1 pairing interaction and surface-peaked multipole- 
multipole interaction terms (quadrupole, octupole and 
hexadecupole) as in Ref. Q. Since ^^Cu is an iV = Z (i.e.. 
To — 0) odd-odd nucleus, its T = 1 levels have analog 
levels in its neighboring even-even nuclei with = ±1. 
The lowest T = and T = 1 levels can be quite close be- 
cause of the pairing energy. The experimentally observed 
ground state of ^*Cu has T = 0, while the lowest T = 1 
state is close to it with an excitation energy of Ex{T = 
1) = 0.203 MeV. However, for the good-sign Hamiltonian 
of Ref. @, we find E{T = 1) to be weU below E{T = 0). 
The energy E{T = 1), describing the ground state of the 
even-even nucleus with T^ = 1, is determined directly 
from {H)a,t^=i = {Tta,t^=i{HU,)/Tva,t^=i{U,))w at 
low temperatures (as in Ref. [9|), while E{T = 0) is de- 
termined from Eq. (O with Tq = and O ^ H. All 
other energies E{T) (for T > 1) are determined from 
{H)a,t^=t, similarly to E{T = 1). Once we determine 
the energies E{T), we adjust their differences to match 
the experimental values of Ex{T). 

In Fig. [T] we present the T- and 7r-projected level den- 
sities of ^^Cu with the corrected values of E{T). An em- 
pirical BBF for the T-projected level densities (including 
the Lang-LeCouteur modification [4]) is given by 



Pt{E,) = g{2T + l)^a'^/\E, - A + tyV^ 



exp 2^a[E^ - A) 



(9) 
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FIG. 1: TTT-projected level densities of ^^Cu, obtained from 
SMMC calculations and shifted to match the experimental 
values of Ex{T). Symbols are as follows: solid (open) circles 
for T — 0, n = + ( — ), solid (open) triangles for T = 1, tt = + 
( — ), solid (open) diamonds for the T — 2, n = + ( — ), and 
solid (open) inverted triangles for T = 3, tt = + ( — ) densities. 
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FIG. 2: Total SMMC level densities of ''^Cu. Solid squares 
represent the density obtained from the sum of the Ttt- 
projected densities after the lowest energy for each isospin 
T is corrected by 5E{T), while the open squares describe the 
total density calculated from unprotected SMMC. Lines are 
fit to the BBF ([T]). Crosses are the results of the perturbative 
aT^ correction with a = 1.2 MeV 



where g = 1 for a T7r-projected density and g — 2 when 
parity is not projected. The projected SMMC level den- 
sities can be well fitted to Eq. ^ if we use T7r-dependent 
values for a and A. The correction to E{T) is described 
by a shift of the corresponding A, but does not affect a. 

The total level density is obtained by summing over the 
TTT-projected level densities (after taking into account the 
shifts 5E{T)). The corresponding total level density for 
^*Cu is shown in Fig. [5] (solid squares). For comparison 
we also show the total level density obtained by a direct 
SMMC calculation (i.e., without T projection) for the in- 
teraction in Ref. [8] (open squares). Fitting the SMMC 
density to Eq. H]), we obtain a = 5.550 ± 0.021 MeV"^ 
and A = -0.658 ±0.076 MeV after projection and SE{T) 
correction, while a = 5.827 ± 0.022 MeV"^ and A 
1.069 ± 0.078 MeV for the unprojected density. We ob- 
serve that the inclusion of the proper isospin dependence 
in the Hamiltonian leads to a significant enhancement of 
the total level density at low excitation energies. This 
can be understood as follows. In the unprojected SMMC 
calculation, the ground state has T = 1 and the level 
density at low energies is dominated by its T = 1 compo- 
nent. The pairing energy in the T — 1 component leads 
to a large A and thus too low level density at low en- 
ergies. Once we include the proper isospin dependence 
(using the T projection method), the ground state has 
T = and the level density at low energies is enhanced. 

We also show in Fig. [2] the total density obtained from 
the perturbative aT^ correction 15. 1^ (crosses). In this 
approach, the energy shifts of the T > To{= 0) compo- 



nents are assumed to be small. We find that the total 
level density enhancement at low excitation energies to 
be too small in the perturbative approach. 

Figure [3] shows the T- and 7r-projected level densities 
of the even-even N = Z nucleus ^°Zn after correcting the 
values of E{T). In Fig. 2] we compare the total density 
obtained from the sum of the projected densities with 
the total level density determined from a direct SMMC 
calculation (with no SE{T) shifts). A fit to Eq. (P) gives 
a = 5.947 ± 0.020 MeV"^ and A = 2.102 ± 0.053 MeV 
for the corrected density (a = 6.185 ± 0.009 MeV~^ and 
A = 1.679 ± 0.059 MeV for the unprojected density). In 
this nucleus the ground state has T = both experimen- 
tally and in the unprojected SMMC calculation. Hence 
the level density at low energies is dominated by its T = 
component either with or without the 6E(T) corrections. 
Nevertheless, we observe that the SE{T) corrections are 
important at higher energies E^ > 10 MeV, where they 
suppress the total level density. To explain this, we note 
that there are two counteracting effects in the isospin 
dependence of the level density: the {2T + 1) factor in 
Eq. ([9]) tends to enhance the projected density for larger 
values of T while the increase of A with T tends to sup- 
press the projected density at larger T. We find that the 
T = 1 component is dominant at E^ > 10 MeV in the 
unprojected total density. In contrast, once we take into 
account the proper energy shift 5E{T = 1) w 1.5 MeV, 
we find the T = 1 density to be substantially suppressed 
and comparable to the T ~ Q density for E^ > 10 MeV, 
as is seen in Fig. [31 This results in an overall suppression 
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SMMC method. Symbols and lines are as in Fig. [T] 



of the total level density for the higher energy regime. 
For this nucleus the perturbative method gives almost 
the same level density as the unprojected density. 

While the isospin-dependent corrections are important 
in = nuclei (and possibly in some jT^j — 1 odd- 
odd nuclei), we have confirmed that these corrections are 
insignificant for \Tz\ > nuclei (at least up to E^: ~ 
20MeV). Since levels with T > Tq (Tq > 1) are well 
separated in energy from the T = Tq levels already with 
the good-sign Hamiltonian, the density of all T > Tq 
levels is less than half the density of the T — Tq levels. 

In summary, we have developed an efficient isospin pro- 
jection method for an isospin-conscrving Hamiltonian in 
the SMMC approach and applied it to nuclei in the com- 
plete pf + 59/2 shell. This isospin projection is exact 
sample by sample and thus leads to statistical errors 
that are typically smaller than the statistical errors in 
a corresponding spin projection method. We have used 
this projection method to take into account the proper 
isospin dependence of the nuclear interaction, avoiding a 
sign problem that occurs when such an isospin-dependent 
interaction is included in unprojected calculations. For 
N = Z nuclei, we find that this isospin dependence can 
lead to significant corrections to the total level density. 

This work is supported in part by the U.S. DOE 
grant No. DE-FG-0291-ER-40608 and as Grant-in-Aid 
for Scientific Research (C), No. 19540262, by the MEXT, 
Japan. Computations were carried out on the PC cluster 
Helios and IBM SP3 in JAERI, and on CP-PACS at the 
Center for Computational Physics at the University of 
Tsukuba. 
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FIG. 4: SMMC total level densities of ""Zn. Symbols are 
as in Fig. [2] except that the perturbative result is not shown 
here. 
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